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A Dyson hierarchical model for Anderson localization, containing non-random hierarchical hop- 
pings and random on-site energies, has been studied in the mathematical literature since its intro- 
duction by Bovier [J. Stat. Phys. 59, 745 (1990)], with the conclusion that this model is always in 
the localized phase. Here we show that if one introduces alternating signs in the hoppings along the 
hierarchy (instead of choosing all hoppings of the same sign), it is possible to reach an Anderson 
localization critical point presenting multifractal eigenfunctions and intermediate spectral statistics. 
The advantage of this model is that one can write exact renormalization equations for some observ- 
^-—^ . ables. In particular, we obtain that the renormalized on-site energies have the Cauchy distributions 

for exact fixed points. Another output of this renormalization analysis is that the typical exponent 
of critical eigenfunctions is always atyp = 2, independently of the disorder strength. We present 
, numerical results concerning the whole multifractal spectrum f{a) and the compressibility x of the 

level statistics, both for the box and the Cauchy distributions of the random on-site energies. We 
discuss the similarities and differences with the ensemble of ultrametric random matrices introduced 
recently by Fyodorov, Ossipov and Rodriguez [J. Stat. Mech. L12001 (2009)]. 
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I. INTRODUCTION 



To better understand the notion of phase transition in statistical physics, Dyson [T| has introduced long ago a 
hierarchical ferromagnetic spin model, which can be studied via exact renormalization for probability distributions. 
Since the hierarchical couplings correspond to long-ranged power-law couplings in real space, phase transitions are 
possible already in one dimension. This type of hierarchical model has thus attracted a great interest in statistical 
I physics, both among mathematicians [2t5| and among physicists [6 8]. In the field of quenched disordered models, 
' hierarchical models have also been introduced for spin systems with random fields or with random couplings 
I p^| - [T2j . as well as for Anderson localizati on jl3l - [T9| . on which we focus in this paper. 
Q ■ In the context of Anderson localization (2C|. the hierarchical models p"3l - [l9j contain long-ranged hoppings decaying 
O ' as a power-law of the distance in real space. To discuss the possibility of an Anderson localization transition, it is thus 
useful to recall first what is known for non- hierarchical long-ranged models. For Anderson models with long-ranged 
' hoppings V{r) presenting the typical asymptotic decay as a function of the distance r 

V 

^ : no « - (1) 

, one expects that Anderson critical points occur in dimension d for a critical value of the exponent given by 

g : a, = d (2) 

whereas a > d correspond to the localized phase (power-law localization) and a < d corresponds to the delocalized 
phase (see the review 21] and references therein). These critical points have been mostly studied in dimension 
d = 1 within the PRBM model (power-law random banded matrices model) i22-.34]. in particular from the point 
of view of their multifractal spectra (studies in d > 1 can be found in (ssMstI ). In hierarchical models, one thus 
also expects the same criticality criterion as in Eq. [51 In the 'ultrametric random matrices ensemble' introduced in 
[l7| , it has been found that criticality indeed corresponds to ac = d = 1 (Eq. [5]) and multifractality properties of 
eigenfunctions have been studied llTHlQf . In particular, in the so-called 'strong multifractality' regime, the powerful 
Levitov renormalization method (Sst (se e also the reformulation as some type of 'virial expansion' in Refs [3914421 ]) can 
be used both for the PRBM model [2^ and for the 'ultrametric random matrices ensemble' [l3| and yield the same 
leading order result, which can also be derived via a simpler direct perturbation theory for eigenstates [43| . However 
the 'ultrametric random matrices ensemble' considered by physicists [l7l - [l9j is not the only possible way to define 
a Dyson hierarchical model for Anderson localization, and another model containing random on-site energies and 
non-random hierarchical hoppings, has actually been considered previously by mathematicians [1343 j some 
region of parameters that does not allow to reach the critical point ac = d = 1. As a consequence, their results concern 
the properties of the localized phase. The aim of the present paper is to show that it is possible to reach the critical 
point (Tc = d = 1 in this tyTpe of models, if one chooses properly some parameter at a well-defined negative value 
(whereas mathematicians [l3l - [la ] seem to have always considered that this parameter was positive). The difference 
between the two types of hierarchical models is as follows : in the model considered in |13l4l6l ] and in the present 
paper, it is directly the matrix of non-random hoppings that presents a hierarchical block structure (see more details 
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in section |TT] below), whereas in References |T7l - [T9| . it is only the matrix of the variances of the random hoppings that 
presents a hierarchical block structure. Note that in the context of spin models with random couplings, these two 
types of hierarchical models have been also introduced, with either a hierarchical structure of the couplings [HI 
or a hierarchical structure of the variances of the random couplings fl^ . 

The paper is organized as follows. In section |TT1 we introduce the 'Dyson hierarchical model with random on-site 
energies and non-random long-ranged hoppings' that we consider in this paper, and we describe an exact rcnormaliza- 
tion procedure in each disordered sample. In section Hill we analyse this renormalization procedure for the pure case, 
we describe the properties of the energies and eigenstates as a function of the hierarchical parameter, and explain 
how to obtain a physical model presenting a real-space power-law exponent CTc = 1. We then turn to the disordered 
model. In section HVl we analyse the renormalization of the on-site energies, we find an exact solution for the Cauchy 
distribution, and we present numerical results for the box disorder. In section [Vj we describe how the two-point 
transmission can be computed via renormalization, and we derive the typical exponent Ktyp = 2 independently of the 
disorder strength. In section IVll we translate the previous result into the typical exponent atyp — 2 for eigenstates, 
and show numerical results for the singularity spectrum f{a) as a function of the disorder strength W. In section IVIll 
we describe our numerical results concerning the intermediate statistics of energy levels. The anomalous weak-disorder 
regime is discussed in section IVIIII Our conclusions are summarized in section HXl 

II. DYSON HIERARCHICAL MODEL FOR ANDERSON LOCALIZATION 

A. Definition of tiie model 




^1 ^2 ^3 ^4 ^5 ^6 ^7 ^8 

FIG. 1: Dyson hierarchical model of Anderson localization defined by the Hamiltonian of Eq. |3]with = 3 and 2^ = 8 sites. 
Each site i is characterized by a random on-site energy e^. Two sites i and j are related by a non-random hopping Hij — V„(ij) 
where n{i, j) represents the generation n where i and j are related in the binary tree hierarchical structure. For instance, the 
site i = 1 is connected to the site i = 2 via the hopping Vi, to the sites i = 3, 4 via the hopping V2, and to the sites i = 5, 6, 7, 8 
via the hopping V3. 

We consider the following Anderson tight-binding model defined for L = 2^ sites by the Hamiltonian (see Fig. [T]) 

i/^[{ei,...,e2«}; {Vi, ..Vn}] ^Y.^^\i >< i\ (3) 

1=1 

+Vi [(|1 >< 2| + |2 >< 1|) + (|3 >< 4| -I- |4 >< 3|) + ...] 

+V2 [(|1 > +|2 >) (< 3|+ < 4|) + h.c. + (|5 > +|6 >) (< 7|+ < 8|) + h.c....] 

-1-^3 [(|1 > +|2 > +|3 > +|4 >) (< 5|+ < 6|+ < 7|+ < 8|) + h.c. + ...] 

+Vn [(|1 > +|2 > +... + |2"-i >) (< 2^-1 + 1|+ < 2^-1 + 2| + ...+ < 2^1) + h.c] 

where the are the independent on-site random energies drawn with some distribution, and where the parameters 
(^1, V2, ...) represent the non-random hoppings at different levels of the hierarchy. In Dyson hierarchical models, it is 
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usual to take couplings that depend exponentially on the generation n 

Vn = V^i7"^' (4) 

To make the link with the physics of long-ranged one-dimensional models, it is convenient to consider that the sites 
i of the Dyson model are displayed on a one-dimensional lattice, with a lattice spacing unity. Then the site i = 1 is 
coupled via the hopping Vn to the sites 2^~^ < i < 2". At the scaling level, the hierarchical model is thus somewhat 
equivalent to the following power-law dependence in the real-space distance L„ = 2" 

|K|= i^^r'^^ with a(7) = -^ (5) 

It seems physically natural to require here that these hoppings decay with the distance L„ with a positive exponent 
cr(7) > corresponding to I7I < 1, i.e. 

- 1< 7 < 1 (6) 

We will see later in this paper that other physical requirements may lead to other restrictions in the choice of the 
parameter 7 (see section IIlID[) . 



B. Exact renormalization procedure 

1. Change of basis 

From the form of the Hamiltonian of Eq. [31 it is clear that it is useful to begin with a change of basis 
from the real space natural basis {|1 >, |2 >, |3 >, |4 >, |2^ - 1 >, |2^ >} 

to the new basis {|1+ >, |1_ >, |2+ >, |2_ >, |2^"^ >, |2^~^ >}, with the following notation for 1 < i < 2^"! 



\i+ > 



|2i - 1 > +\2i > 



In this new basis, the Hamiltonian of Eq. [3] reads 

iJjv[{ei,...,e2«}; {Vi,..Vn}]^ J2 >< i+\ + >< i-\ + V^+_ >< i-\ + \i- >< (8) 

i=l 

+2V2 [|1+ >< 2+1 + |2+ >< 1+1 + |3+ >< 4+1 + |4+ >< 3+1 + ...] 

<A+\) + h.c. + ...] 

N-2^^f^roN-2 , , o^ , , , <2^-'\)+h.c.] 

in terms of the new parameters 

(^^^ 1 -I- 

+ Vi 

V.,. (9) 





- h.c. + 


>)( 


;< (2^- 




-1 + 




2 


£24- 


-1 + ^2i 




2 


e2i~ 


1 — £21 



2. Elementary renormalization step 

For Anderson localization models, Aoki [43 - l46| has introduced an exact real-space renormalization procedure at 
fixed energy which preserves the Green functions of the remaining sites. This procedure has been further studied 
for one-particle models in [47l - l49| . It has been extended in configuration space for two-particle models [s^l and for 
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manybody localization models [5lj. It can be also used in other physical contexts, like phonons in random elastic 
networks ^52i] . Let us now explain how it works for our present purposes. 

From the expression of the Hamiltonian in the new basis of Eq. [HI it is clear that for each 1 < i < 2^^^, the vector 
\i- > is only coupled to its symmetric partner |i+ >. As a consequence, if one considers the Schrodinger equation at 
some given energy E 

E\tP>^HN^> (10) 
the projection onto < with the usual notation < i±\ip >= V'(*±) 

Ei^ii^) - e,_7^(z_) + V,+-ijii+) (11) 
can be used to eliminate ipi}-) in terms of tp{i^). One obtains for > the renormalized on-site energy 



new / 77i\ _ , ' iH — 

ei+ (^) - ei+ + 



£21-1 + £21 

2 



f e2t-l-e2i \ 

V 2 ; 



+ Vi + ^ r (12) 



E- 

For the remaining degrees of freedom {|1+ >, |2-|_ >, |3+ >, |2^^^ >}, the Hamiltonian reads 

2"-i 

i?^-™[{ei,...,e2«};{v^i,..vv}] = E errK+><* + i 

n=l 

+2V2 [(|1+ >< 2+1 + /i.e.) + (|3+ >< 4+1 + h.c) ....] 
+2V3 [(|1+ > +|2+ >) (< 3+1+ < 4+1) + h.c. + ...] 

+2Vn [(|1+ > +|2+ > +... + \2^-^ >) (< (2^-2 + 1)^1 + < 2^-i|) + h.. 
- i/Ar-i[{e?r, e^^"^!}; {V( = 2V2, .., = 2Vn}] (13) 

In summary, the Hamiltonian defined for the 2^ initial sites with the random energies (ei, £2") ^'iid the 
hierarchical hoppings (Vi, V2, ...Vat) leads, after the decimation of the 2^^^ antisymmetric combinations {|1_ > 
,|2_ >,..., |2^-i >},to the same form of Hamiltonian for the remaining 2 symmetric combinations {|1+ >, |2+ > 
, |2+ >}, with renormalized parameters : the renormalized energies are given by the rule of Eq. 1121 whereas the 
hierarchical couplings are simply given by V{ — 2V2, .., = 2Vn. 

Before we consider the disordered case, let us first describe in the next section the properties of the pure model. 



.c. 



III. PROPERTIES OF THE PURE MODEL = eo 



In this section, we describe the properties of the pure model, where all on-site energies in the Hamiltonian of Eq. 
[3] take the same value = cq, and where the hierarchical couplings follow the geometric form of Eq. |4l 

A. Enumeration of eigenstates 

We follow the renormalization procedure introduced in the previous section fllB I In the pure case, the parameters 
after the first change of basis do not depend on i and simply read (Eq. [S]) 

=eo + Vi 

Vl"-} =0 (14) 

Since the intercoupling vanishes V_l^2 = 0, the 2^^^ states \i^ > are completely decoupled : they are exact eigenstates 

of energies e^''. So in the pure case, the decimation step is completely trivial, and there is no renormalization of the 
on-site energies for the remaining states |i+ >, ( Eq. [T7\ reduces to the identity e^^^iE) = e^^). 
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The second renormalization step is thus similar, but with the new values (using V( = 2V2 from Eq. I13p 

ef - e^l^ + Vi = eo + 14 + 2^2 



Vi = eo + V^i - 2V^2 (15) 



So by iteration, the n-th renormalization step is again similar, but with the new values 

ei"^ = eV""'^ + 2"-iK„ = 60 + 1^1 + 2V^2 + ... + 2"-V„_i + 2"^ V„ 

eL"^ = eV""'^ - 2"-iK = eo + V^i + 2^2 + ... + 2"-V„_i - 2"-!^ (16) 
When the hierarchical couplings Vn follow the geometric form of Eq. 21 these energies simply read 



eV"' =eo + yi 5^(27)'-^ = £0 + ^1^ 



27 

» _>-i) on-iT. , l-(27)"-i2(l-7) 

1 - 27 



1=1 

\ n — 1 r 



e'.' ^ eV"^^ - 2"- V„ ^ 60 + ^1 ^ V oJ ^ (17) 



We may now enumerate all the eigenstates, or more precisely all the subspaces associated to eigenvalues. 

1. Subspace of dimension 2^~^ associated with the energy 

We have seen above that the subspace associated to the energy can be constructed from the basis (Eq[7]) 

with i = 1,2, 2^^^. Each vector |z_ >(^) is thus extremely 'localized' : its support is made of two sites only. This 
property may seem very strange by comparison with non-hierarchical translation invariant models, where eigenstates 
are given by plane-waves that are completely delocalized over the whole lattice. However here also we can make 
plane- wave linear combinations of the basis of Eq. [TS] to obtain a new-basis (noting j = V— 1) 

2Af-l 

l'?>^'^=7^E^^'"N->^'^ (19) 

with the wave-vector q — m{2n/2^~^) where m takes the integer values m — 0,1, .., 2^~^ — 1. Each vector \q >(^) is 
then delocalized on the whole lattice, with a constant weight | < i\q >^^^ ^ = 1/2^. 

2. Subspace of dimension 2^~^ associated with the energy e^' 

(2) 

The subspace associated to the energy e_' can be constructed from the basis 

^(2)^ |(2»-l)+>(i)-|(2»)+>(i) ^ |4z-3>+|4»-2>-|4z-l>~|4»> 
' V2 2 ^ ' 

with i = 1,2, ...,2^~^. Each vector |i_ > is 'localized' on four sites only. But again, one may use instead a basis of 
delocalized plane-waves 



^ y e^"«'|»_ >(2) (21) 
V2^ 



with the wave-vector q = m(27r/2^ ^) where m takes the integer values m = 0, 1, .., 2^ ^ — 1. 

It is now clear how this construction can be pursued, and we only mention the last two generations. 
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3. Subspace of dimension 2 associated with the energy 
The subspace associated to the energy e^^^-* can be constructed from the basis 

II- 'fg I. >^ 

/nN-1 I oiV-2 

|o U ^ 7V-2 ^ /2 +2 ^ 



iV-1 



(22) 



Subspaces of dimension 1 associated with tL^ and 

The last renormahzation step yields that the energies ef^-* and e^i^"* are non-degenerate and are associated respec- 
tively with the eigenstates 

'2™-^ 2 

and 

2 

72 - 72^' 

One may check that the total number of states is 2^ as it should 



i=l z=2"-i + l 



2^ - 1 

[2^-1 + 2^-2 + ... + 2 + 1] + 1 = ^— ^ + 1 = 2^ (25) 



B. Green function of the pure case 



From this enumeration of eigenstates, one obtains directly that the Green function at coinciding points reads 

12 ^ 1 1 



,20, 

a " n=l Z e_ ^ 

(because at the given point 1, only one state of each generation gives a non-zero weight |'0ct(l)p). 

More generally, the Green function between two distinct points (1,2^) has a similar form as Eq. except that 
the sum begins only with the state of generation p which is non- vanishing both at 1 and 2'' 

^ , ^_V V'a(l)C(2^) _ V ^ , W (OT, 

- 2. ,_E^ - 737^ + JTTlv) (27) 

a n=p ^ t_ ^ 

C. Density of states of the pure case 

From the Green function at coinciding points of Eq. 1261 one obtains that the density of states in the thermodynamic 
limit — > +00 reads 
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(the normalization is J dEp[E) = X^n^^ ^ = 1 as it should). So in contrast to usual non-hierarchical models in finite 
dimensions that are characterized by a continuum of non-degenerate delocalized states, the density of states of the 
hierarchical model in the thermodynamic limit remains a sum of delta functions, with highly degenerate energies. 
Since the levels are given by Eq. [T7] 



f}:' = €0 + I'l [1 + 27 + ^- + (27)"-' - (27)"-'] 

^.o + >^. '-'Y:;'"~"' PS) 

various behaviors are possible as a function of the hierarchical parameter 7 introduced in Eq. |4] (recall the condition 
I7I < 1 of Eq. El): 

(i) for 1/2 < I7I < 1 : the absolute values of the energies |e!^''| grow exponentially in n. This case is thus unphysical, 
since higher levels of the hierarchy that corresponds to more extended states are associated with higher and higher 
energies. 

(ii) for I7I < 1/2 : as n grows, the energies accumulate near the finite accumulation point 

e!.- = eo + V^iY^ (30) 

(iii) for the special case 7=1/2: the absolute values of the energies le^''] grow linearly in n 



e 



(«) 



eo + Vi[n- 2] (31) 



This case is thus unphysical as the case (i) discussed above. 

(iv) for the spec 
in contrast to (iii) 



(iv) for the special case 7 = — 1/2 : the energies e^"* take alternatively two values, and the model is thus 'physical' 



D. Physical region of the model 

From the above discussion, the requirement to have a bounded density of states for the discrete model of Eq. [3] 
with the hierarchical couplings of Eq. |4] yields the following domain for the hierarchical parameter 7 

4<7<^ (32) 

instead of Eq. [6l 

From the point of view of the real-space power-law exponent cr(7) of Eq. [Sj our conclusion is thus that the domain 
(7(7) < 1 corresponding to 1/2 < I7I < 1 (case (i) above) is not physical, whereas the domain (7(7) > 1 corresponding 
to I7I < 1/2 (case (ii) above) is physical. Finally the critical value ctc = 1 of Eq. [2] cannot be made physical with the 
choice 7=1/2 (case (iii) above), but can be made physical with the choice 7 = —1/2 (case (iv) above). 

In the presence of disorder, the criterion of criticality for Anderson localization in the presence of long-ranged 
hoppings that we have recalled in the introduction around Eq. [5] leads to the following conclusion : all cases with 
— ^ < 7 < ^ are expected to be in the localized phase, whereas the case 

7c = (33) 

is expected to be critical. In the remaining of this paper, we thus focus on this case, which has not been considered 
previously in the mathematical literature as we now explain. 



E. Relations with the models considered in the mathematical literature 

1. Relation with the model of Ref fldll 

Eq. [26]is in agreement with formula (2.25) of Bovier [13] in the thermodynamic limit N -> -f 00, with a simple shift 
n — >■ n -I- 1 in the labelling of e*^'"* given in Eq. [iTl with the following correspondence of notations : Bovier's model 
corresponds to the parameters eq = 0, 14 = — 1/2 and 



a = 47 



(34) 
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Bovier has assumed a > 0, in order to interpret the pure Hamihonian Hq = — Aq, as a 'hierarchical Laplacian'. 
However if one defines the model via a tight-binding quantum Hamiltonian as we did in Eq. [31 there is no reason 
from a physical point of view to impose a priori some sign constraints on the hoppings T^, i.e. on a = Aj. So the 
interesting critical case of Eq. [551 corresponds to a = —2 in the notations of (Tsj . 

2. Relation with the model of Refs fTMldl l 

In Refs p^ - [l6j , the notion of 'hierarchical Laplacian' has been defined in terms of a sequence of positive numbers 
p„ > for n = 1, 2, .. such that 

-t-oo 

^p„ = l (35) 

n=l 

Translated into our present notations with the quantum Hamiltonian of Eq. [31 this version of the model corresponds 
to the pure on-site energy 

-t-oo 

^0 = E (36) 

n=l 



and to the hierarchical hoppings 



so that the eigenenergies of Eq [THl reads 



+ 00 

^" = E (37) 

m—n 



eL"='^ = £0 - ^1 = 
eL"='^ =eo + Vi- 2V2 = pi 

n-l 

eL"^ =eo + V^i+21/2 + ... + 2"-2y„_^_2"-V„= Ep- (38) 

m— 1 



where the accumulation point is fixed to (Eq. I35p 



e°^^ = E^'" = l (39) 



71=1 

Here again, Refs [IM^ 

impose the positivity of the numbers Pn > to have an interpretation in terms of some 
classical random walk for the hierarchical Laplacian, whereas if one defines the model via a tight-binding quantum 
Hamiltonian as we did in Eq. [3l there is no reason to impose the positivity of all the p„, and the only physical 
constraint is on the eigenvalues e^' that should remain bounded, i.e. the constraint is that all partial sums X]m^\ Pm 
should remain bounded. In particular, the interesting critical case of Eq. [33l corresponds to the choice p„ = (— 1)"^^. 

F. Special case 7c = —1/2 corresponding to a real-space power-law of exponent CTc — 1 

For the special case jc — —1/2, there exist an exact invariance after two steps of RG, and Eg [TTI become 

^> + ,40) 
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As a consequence, in the enumeration of eigenstates of section UlI Al the states (— ) have only two possible energies, 
corresponding to odd and even generations 



whereas the last state (+) has for energy 



^even ^ ^ ^Vl (41) 



ef ^"^"^ = £0 (42) 



In the thermodynamic hmit — > +C50, the Green function of Eg 1271 becomes 



G, 2P (z) = y = + - — — (43) 

n—p ^ ^— ^ ^ 



with the coefficients 



oddf^p. d{2k + l>p) 
9 (2*^) = ^2fcTI 



k=0 



/c=2 

that decay as 1/2^ = 1/Lp with respect to the distance Lp = TP . 

In particular at coinciding points corresponding to p = [2^ = l)i one obtains the simple values 



^^''^"(1) -E^:-^ (45) 
so that the density of states of Eq. [^5] is simply given by the following sum of two delta peaks 

PpureiE) = y{E- t"^") + i(5 (S - e™™) (46) 

This pure model is thus extremely degenerate. As a consequence, the effect of a 'weak disorder' is expected to be 
anomalous at the perturbative level (see more details in section IVIIip . 

IV. ANALYSIS OF THE RENORMALIZATION OF THE ON-SITE ENERGIES 

In this section, we analyse the renormalization procedure of section lll BI when the initial on-site energies are random 
variables. We first consider the case of the Cauchy distribution that leads to an explicit exact solution. We then 
study numerically the renormalization flow starting from a box distribution. 

A. Cauchy disorder : exact solution 

1. Reminder on some specific properties of the Cauchy distribution 

As is well-known, the Cauchy distribution is the only probability distribution which is stable both by addition and 
by inversion. More precisely, if one denotes the Cauchy distribution of mean a and width b by 

CaA^) ^ ^ .1 (47) 
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one has the two properties : 

(i) addition : if xi is distributed with Ca-i^biixi) and X2 is distributed with Ca2.b2i^2)i then x — X1+X2 is distributed 
with Ca,bix) of parameters a = ai + 02 and 6 = 61 + 62- 

(ii) inversion : if a; is distributed with Ca.tix), then y = 1/a; is distributed with CA.siy) of parameters A = a/(a'^+b'^) 
and 5 = 6/(0^ + 6^)- 



2. Renormalization flow of the on-site energies ei± 

The RG Equation [T2I that describes how the on-site energies evolve upon the first RG step (Eq. [T3|) . can be 
rewritten as 

1 _ 1 1 

Let us recall the notations used in this relation : the energy E and and the hopping Vi are parameters, €2^-1 and £2; 
are independent on-site energies of the initial model, and 6^+ is a renormalized on-site energy after one renormalization 
step. 

We now apply the properties of the Cauchy distribution recalled above in section HV A 11 We assume that the two 
independent on-site energies 62^-1 and £2?: of the first generation are drawn with the distribution Ca^.tii^) with the 
values 

oi = eo 

bi^W (49) 

where W represents the disorder strength around the averaged value eo- 
As a consequence, one obtains that 

(1) (e2i-i — E — Vi) and (e2i — E — Vi) are distributed with Cai-E-Vi.bi 

(2) l/(e2i-i — E — Vi) and l/(e2i — E — Vi) are distributed with Ca,b of parameters 



{ai^E~Viy + bi 



(3) the sum S = l/(e2j:-i - E - Vi) + l/{e2i - E -Vi) is distributed with C'2A,2B- 

(4) the inverse / = l/S* is distributed with Caj,bi of parameters 

2A ai-E -Vi 

ai 



(2A)2 + (2B)2 2 

" (2A)2 + {2By ' y ^^^^ 

(5) the on-site energy e"_J:™ ^ E + 2Vi + 21 after one RG step is then distributed with Ca2,b2 of parameters 

a2 = E + 2V1+ 2ai =01 + ^1 

62 = 2bi = 61 (52) 

i.e. the width of the distribution is conserved, whereas the averaged value in shifted by Vi. 

Taking into account the evolution of the lowest coupling after (n — 1) renormalization steps (see Eq. I13|) we obtain 
that if the on-site energies of the initial model are distributed with the Cauchy distribution C'ai^bi with Eq. 1491 then 
the on-site energies after n renormalization steps are distributed with the Cauchy distribution C'a„.b„ of parameters 

a„+i = a„ + V}''-'^ =eo + + V^'K.. + v/""'^ =01 + ^1+ 2^2 + 2V3 + ... + 2"-^^ 

bn+i =b,, = bi^W (53) 

For W = 0, we recover of course the pure model analyzed previously in section Hill In particular, the physical cases 
where the on-site energies do not flow to infinity correspond to the domain of Eq. 1321 
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3. Critical case fc = —1/2 ; exact invariance after two RG steps 

As previously discussed in section UlI Fl concerning the pure case, the special case 7c — —1/2 corresponding to a 
real-space power-law of exponent cr^ = 1 has the special property to be exactly invariant after two RG steps. This 
property survives in the disordered case for the Cauchy distribution : the parameters of the Cauchy distribution of 
Eq. [53] are stable after two renormalization steps 

a2n+l = 0-odd = 0-1 ~ f-O 

a2n = eleven = Oi + Vi = + Vi (54) 



B. Case of other initial on-site energy distributions at criticality "fc = —1/2 




FIG. 2: Integrated distribution In{E) = deP„{e) of the probability distribution -Pn(e) of the renormalized on-site energies 



obtained via the RG rule of Eq. [56] as a function of the generation n = 2, 4, 6, 8, 10 starting at n = from the box distribution 

■ fc ^ 

' ■K[(x-a)^+b^\ 



of Eq. [55] with VF = 5 : the convergence towards the integrated Cauchy distribution Icauchy{E) — J^_^dx- * 



^ -I- ^ arctan (^ ^ ^ ° ) is very fast. 

As explained in the introduction around Eq. [21 one expects that a model with real-space power-law hopping a = 1 
will be critical for any on-site energy distribution. In our present case where tic = 1 is realized for the value 7c = —1/2, 
we thus expect criticality to occur for any initial distribution of the on-site energies. It is then interesting to study 
the renormalization flow of the on-site energies. As an example, we have considered the case where the initial on-site 
energies are drawn with the box distribution 

and we have computed numerically the probability distribution of renormalized on-site energies obtained by the RG 
equation for 7c = —1/2 

1 - 1 I 1 (56) 

M::^„y^(„i)n-i (4ri^-i?)-yi(-i)«-i (4r'^-i?)-vi(-i)"-i 

As an example, we show on Fig. [2] the numerical results corresponding Vi = 1, E = starting from the initial 
distribution of Eq. [55] with IF = 10 : after a small number n = 2,4,6,8, 10 of RG steps, we obtain a very rapid 
convergence towards a Cauchy-like distribution. Our conclusion is thus that the Cauchy distributions found above 
seem to be the only attractive fixed points. 



V. RENORMALIZATION PROCEDURE TO STUDY THE TRANSMISSION PROPERTIES 



In quantum coherent problems, a convenient characterization of transport properties consists in defining a scattering 
problem where the disordered sample is linked to incoming wires and outgoing wires, and in studying the refiection and 
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transmission coefBcients. This scatterin g th eory definition of transport, first introduced by Landauer [56{ . has been 
often used for one-dimensional systems [57H59| and has been generalized to higher dimensionalities and multi-probe 
measurements (see the review [60|). 

In this section, we assume that the first site i = 1 is connected to an incoming wire, and that the last site i — 2^ 
is connected to an outgoing wire (see Figure [Sj . 



Elementary RG step for the couplings to the external wires 



We call A the point of the wire of on-site energy ca that is connected to the site i = 1 via the hopping Va,i- During 
the change of basis of Eq. [71 the hopping between A and 1 becomes 



Va {\A >< 1| + |1 >< A\) = Va+ {\A >< 1+1 + |1+ >< A\) + Va- {\A >< 1-1 + |1_ >< A\) 



with 



Va+ - Va- 



Va 
V2 



Since |1_ > is only coupled to |1-|_ > and \A >, the Schrodinger equation at energy E 



(57) 



(58) 



(59) 



can be used to eliminate tp{i-) in terms of and '^{A), so that in addition to Eqs[T2l the renormalized on-site 

energy of site A and the renormalized coupling F^^™ read using Eqs [5] 



^A 



xrnew 



(■A 



Va^ 



E-ei_ 



£A 



2{E+Vl)- (6^+62) 



Va-V-+ ^ Va 
E-ei_ 



ei - £2 



2iE + Vi) - (ei +£2) 



(60) 



In the following, the important point is that the the renormalization rule for the renormalized coupling to A can be 
rewritten as 



l/new 



1 



Va V2 [ 



e2~E^Vi 



Similarly, the hopping Vb between B and the last site 2^ is renormalized during the first RG step into 

e(2"-i) ~ E ~ Vi 



jrnew 

Vb 



1 

V2 



<^(2«-l)+<^2« 



E-Vi 



(61) 



(62) 



13 



B. Effective transmission for a finite system of size Ln = 2^ 

We consider a finite system of size Ljv — 2^. During tire first {N — 1) RG steps, tire hoppiirgs Va and Vb to 
the external wires are renormalized independently. Denoting by £^"2 the first two renormalized on-site energies after 



{n— 1) RG steps (so that the initial on-site energies correspond to ej" — ei and £3" 
renormalized renormalized coupling, one obtains the product structure 



V. 



(JV-l) 



A 



Rn-1---R2Ri 



where 



Rv = —7= 



1 



Ap) 



E + Vi{-lY 



E + Vi{-iy 



£2), and V"^"'' the successive 
(63) 

(64) 



Similarly, denoting by Cj"'' and e^"'' the two last renormalized on-site energies after (n — 1) RG steps (so that the 



initial on-site energies correspond to 



, ~(n=l) 

e2N_i and e-j^ 



V, 



(N-l) 



B 



£2"), the renormalized coupling to site B reads 
RN-1-R2R1 (65) 



where 



Rv 



1 

71 



'^^-E + Vi{-iy 



(66) 



After these (TV— 1) RG steps, the renormalized system contains only two renormalized sites |l(jv-i) > and |2(^_i-) >, 
and the corresponding renormalized Hamiltonian reads 



Mast 



(N-l) 



|l(Af-l) >< 1(7V-1)I + ^^|2(JV-1) >< 2(Ar_i)| 

-K(_l)^-ll/l (|2(^_i) X l(^r_i)| + |l(Ar_i) X 2(^r_i)|) 



-V. 



{N-l) 



^ (|A>< l(jv_i)| + |l(^r_i) >< A|) +1/^ 



{N-l) 



(|S><2(^r_i)| + |2(jv-i) ><B\) 



(67) 



The successive elimination of |l(jv-i) > and |2(-jv-i) > via Aoki RG rules finally yields the following renormalized 
coupling between the two wires ends A and B 



vrj = (-1) 



V{N-1).AN-1) 



{e-A--^\e. 



(68) 



Since the renormalized on-site energies and £ 



(N-l) 



have a fixed probability distribution independent of iV, 



the statistical properties of the Landauer transmission for large TV can be obtained from the the renormalized hoppings 
alone 



Tn - \V, 



e-ffl-2 



2|T.(W-1)|2 



AB 



W, 



{N-l),2 



B 



(69) 



In conclusion, the only important A^-dependent factors are the two independent renormalized couplings 



(N-l) 



and 



V, 



{N-l) 



B 



given by the products of Eqs [SS] and [ 



Multifractal statistics of the transmission 



At Anderson localization transitions, the two-point transmission displays multifractal properties in direct corre- 
spondence with the multifractality of critical eigenstates [H, [H, [53| : the critical probability distribution of Tl over 
the disordered samples takes the form 



Pioh (Tl ^ L-'') dT cx L'^^'^Uk 



(70) 
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and its moments involve non-trivial exponents X{q) 



rpq 



L— >-oo 



(71) 



As stressed in (53| , the physical bound < 1 on the transmission implies that the multifractal spectrum exists only 
for K > 0, and this termination at k = leads to a complete freezing of the moments exponents 



X{q) = X{qsat) for q > qsat 

at the value qsat where the saddle-point of the integral of Eq. [711 vanishes K{q > qsat) = 0. 
In our present case, the moments of the transmission of Eq. 1691 read 



rpq 



2\q 



w. 



A 



2q 



1 



-X{q) 



^ 2-NX(q) 



Using Eq. 



one finally obtains that the exponents X{q) read in terms of the variables Rp of Eq. 

2 



X(q)^- lim 



In 



|i?JV-l...i?2i?lP9 



(72) 



(73) 



(74) 



The exponents X{q) thus represent the 'generalized Lyapunov exponents' for the moments of products of the variables 
i?p, that are correlated via the renormalization procedure. Whereas for arbitrary g, these correlations make difficult 
the analytical computation of X{q), the problem becomes simple if one considers the typical behavior near g = 
where one has to evaluate sums instead of products, as we now describe. 
The typical transmission 



rptyp 



is expected to decay at criticality with some power-law 



typ 



(75) 



(76) 



^ L-¥00 L^typ 

where Ktyp is the value where $(k) is maximum and vanishes ^(Ktyp) = 0. In our present case, we may rewrite Eq. 
1551 as 



In IK 



(™)|2 _ 



A 



(77) 



Using the stationary measures Ca,b with aeven ~ + Vi and aodd — eo for renormalized energies at even and odd 
steps (Eq. IM)) . one obtains that (ei + €2)/ 2 is distributed with the same law as ei and £2- As a consequence, the 
average over the stationary measure yields 



« Ini?: >>= -ln2 



and we obtain the typical decay 



In 1^1") P 



-TO In 2 



In L„ 



in terms of the length L,„ = 2™. Finally, the full transmission of Eq. [SSlhas for typical behavior in L 
h^~ln|yj^"^^|2+ln|\4^"^^|2 =-2(A^-l)ln2--Kf3,plnLAr with Ktyp ^ 2 



(78) 



(79) 



(80) 



Note that this corresponds to the naive estimate T/v ^ ^1/ Lj^ in terms of the direct coupling Vjv between the 
sites 1 and 2^. In the next section, we explain the consequences for the multifractal statistics of eigenstates. 
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VI. MULTIFRACTAL STATISTICS OF EIGENSTATES AT CRITIC ALITY 7^ = -1/2 

One of the most important property of Anderson localization transitions is that critical eigenfunctions are described 
by a singularity spectrum /(a) defined as follows (for more details see for instance the review [2l|V in a sample of 
size L'', the number of points r where the weight |V'(^P scales as L~" behaves as 

UUa) cx (81) 

L— >oo 

The inverse participation ratios (I.P.R.s) can be then rewritten as an integral over a 

Yg{L)= f d'^r\ij{r)\^'^ f da L^^"^ L-'^'' ~ L^^^") (82) 

J^d J L^oo 

where r(g) can be obtained via the Legendre transform formula 

— r(q) = maxo, [f{a) — qa] (83) 



A. Strong disorder regime 

As recalled in the introduction, the 'strong multifractality' regime has been first studied via the powerful Levitov 
renormalization method [s^ (see also the reformulation as some type of 'virial expansion' in Refs [3^-[i3) ■ 
obtained leading order result is the same for the PRBM model [23*1 and for the 'ultrametric random matrices ensemble' 
[TtI [isj , and can also be derived via a simpler direct first-order perturbation theory for eigenstates [i^l • For our present 
model, we may directly apply the analysis of [i^ , where the only important factor is the two-point hopping (see section 
2.1.2. in [4^). We thus conclude that our present model is described by the same 'universal' strong multifractality 
regime for large W 

f^'^^{a) = 1 for < a < 2 (84) 

or equivalently for the exponents T{q) 

T'5J^^(g) =d(2q-l) for g< i 

T^*^(g) =0 for g > i (85) 

The first correction in 1/W discussed in [13, 0; HI; 113 involves an explicit perturbative universal correction for T{q) 
in the region q > 1/2, that should hold also in our present model. But then to obtain the perturbative correction for 
T{q) for q < 1/2, and in particular the behavior of /(a) near the typical value {atyp = a{q — >■ 0)), one has to rely on 
the 'multifractal symmetry' (see more details in section FVl C Sp . Since it is not clear to us whether this symmetry is 
satisfied in the present Dyson model, we cannot conclude that the whole first order correction for all values {a,q) is 
the same as in the PRBM model (23j and in the 'ultrametric random matrices ensemble' [TtI. ITsj. In particular, in 
these two models, the typical value atyp is known to move from afy^ — 2 ai first order (i.e. 2 — atyp oc 1/W), whereas 
in the present model, it seems to us that the typical value atyp remains frozen to the value Oif^jj ~ 2 independently 
of the disorder W, as we now discuss. 



B. Relation with the multifractal statistics of the transmission 



The multifractal spectrum $(k) of Eq. [70] is expected to be directly related to the singularity spectrum /(a) of Eq. 
[Silvia |48|,^,^J 



$(k > 0) = 2 ^f{a ^d+-)-d 
In particular, the typical exponents Ktyp and atyp where ^{Ktyp) = and f{atyp) = d are related by [48l. Issl. [54 



(86) 



atyp = d + 



(87) 
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In our present case where d = 1, the result Ktyp = 2 of Eq. [80l yields 

atyp = 2 (88) 

independently of the disorder strength W. This behavior is very anomalous with respect to other models where 
the typical exponent usually moves continuously between the weak-disorder value atyp{W 0) — d and the strong 
disorder value atyp(W — oo) = 2d. 

In terms of the moments exponents, Eq. [86] is equivalent in the region q < qsat to [H, [s^, [131 

X{q) = 2 [r(g) - d{q - 1)] (89) 



C. Numerical results for the multifractal spectrum as a function of W 





FIG. 5: Cauchy disorder : multifractal statistics of eigenfunctions for the disorder strengths W = 1,5, 10, 20, 40 (a) singularity 
spectrum /(a) (b) corresponding a{q) 
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We have studied via exact diagonalization disordered samples of sizes Ln = 2^ (with N = 7,8, 9, 10, 11, 12 gener- 
ations) with the fohowing corresponding numbers ng^Lj^) of disordered samples 

N = 128; 256; 512; 1024; 2048; 4096 
ns{N) = 51.10^; 11.10^; 215.10^; 326.10^; 3.10^; 780 (90) 

In each disordered sample, we have analyzed the fraction 1/8 of the eigenstates at the center of the spectrum (we 
have checked that the density of states is nearly constant in this region and that the corresponding eigenstates have 
the same statistics of I.P.R.). The multifractal spectrum /(a) is then obtained parametrically in q via the standard 
method of Ref [5^ (see more details in Appendix B of [GlJ). We have chosen the values eo = and Vi = 1. 



1. Results for the box distribution 

We show on Fig. [¥] our numerical results concerning the singularity spectrum /(a) for the disorder strengths 
W = 5,10,20,40,80. At strong disorder W — 80, the multifractal spectrum is very close to the strong multifractal 
universal result of Eq. [84] as expected. As the disorder strength W becomes smaller, the multifractal spectrum 
becomes more curved but keeps a 'strong multifractal' character, with a typical value remaining at the maximal value 
atyp — 2, and an almost vanishing minimal value amin — 0. The inhomogeneity of eigenfunctions is thus always very 
strong. 

For the box disorder below W = 'b, the density of states tend to break into two bands around the two pure delta 
peaks of Eq. [35] that are separated here by ej„g„ — e^^^j = 3Vi = 3. As a consequence, the choice to work around 
_B = is not appropriate anymore at weak disorder, and we discuss this limit in section rVIIII 



2. Results for the Cauchy distribution 

We show on Fig. [5] our numerical results concerning the singularity spectrum f{a) for the disorder strengths 
= 1, 5, 10, 20, 40 : the results are qualitatively similar to the box disorder case. 

Note that for the Cauchy disorder, the density of states can be computed exactly from the pure density of states 

M 

^(^)-/-^^'^^"^^(^ \[(i^-^)^ + i^^] 

i.e. in our present case with the form of Eq. [46] for pP"^'^{E'), we obtain 

2 W I W 

^^^^ ^37r[{E- e^d'^Y + W^] ^ 3 tt [{E - e!!'^-")2 + W^] ^^^^ 

with the numerical values e°^'^ = -1 and e™"" = 2 (Eqs[4T]with eq = and Vi = 1) 



3. Test of the symmetry f{2d — a) = f{a) + d — a 

For any Anderson transition in the so-called 'conventional symmetry classes' [2l|, Mirlin, Fyodorov, Mildenberger 
and Evers [6^ have proposed that the singularity spectrum /(a) of critical eigenfunctions satisfies the remarkable 
exact symmetry 

fi2d - a) = /(a) + d- a (93) 

that relates the regions a < d and a> d. Further discussions can be found in [2l|, [gj, [fiHl ■ In terms of the q- variable, 
the symmetry with respect to the value as — d becomes a symmetry with respect to the value Qs — 1/2 |63| , so that 
one should observe the following fixed point 

"(9-0=1 (94) 

On Fig. [3] (b) for the box disorder case, and on Fig. [5] for the Cauchy disorder case, we find that the curves aw{q) 
for various disorder strength W cross near the point of Eq. 1941 To test more directly the symmetry of Eq. [^3] with 
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FIG. 6: Test of the symmetry of Eg. 1931 : we compare the singularity spectrum /(q) and the function /g (a) = /(2 — a) — (1 — a) 
(a) for the box disorder with W = 5 (b) for the Cauchy disorder with W — 1 



d = 1, we have plotted f{a) and = /(2 — a) — (1 — a) together on Fig. [6l for the box disorder at W — 5 

and for the Cauchy disorder at W = 1 : the difference between the two remains within our numerical errors. For 
larger W, the difference become smaller, as could be expected since the strong multifractality limit of Eq. [Ml satisfies 
the symmetry exactly. In summary, since the deviations with respect to 'strong multifractal limit' are small, the 
deviations from the symmetry are also small, and it seems difficult to obtain a clear numerical conclusion. On the 
other hand, the discussion after Eq. |SS] suggests that the multifractal symmetry is not compatible with our statement 
of Eq. [88] concerning the fixed value of the typical exponent atyp — 2 independently of the disorder strength W. The 
clarification of this point goes beyond the present work. 



VII. COMPRESSIBILITY OF ENERGY LEVELS AT CRITICALITY 7^ 
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FIG. 7: Box disorder : statistics of the number Ml of eigenvalues within a fixed interval /ab = [— A£'/2, +Ai?/2] as a function 
of the system size L. Both the averaged number and the variance grow linearly in L, the ratio of the two slopes yields the 
compressibility x (p) Case W = ^ : we measure x — 0.45 (b) Case = 80 : we measure x — 0.87 

Another important property of Anderson localization transitions is that the statistics of eigenvalues is neither 
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'Poisson' (as in the localized phase) nor 'Random Matrix' (as in the delocalized phase) but 'intermediate' (see for 
instance |66l468l | and references therein ). A convenient parameter is the level compressibility x, which is found to 
satisfy < x < 1 at Anderson transitions (see for instance [l^ and references therein), whereas delocalized states are 
characterized by Xdeioc = and localized states by xioc = 1- 

We have studied via exact diagonalization disordered samples of sizes Ln — 2^ (with = 7, 8, 9, 10, 11, 12 gener- 
ations) with the following corresponding numbers ns{LN) of disordered samples 

Ln = 128; 256; 512; 1024; 2048; 4096 
n,{LN) = 41. 10^; 915. 10"*; 17.10^19. 10"*; 12. 10^; 48. 10^ (95) 

In each disordered sample, we have analyzed the number of eigenvalues within a fixed interval /a_e = +AE/2] 
where the density of states is nearly constant. The averaged number of eigenstates in this interval for disordered 
samples of size L scales as 

< TVl > ~ Lp{0)AE (96) 
The compressibility is then defined by the ratio between the variance and the averaged number 



X (97) 



A. Results for the box distribution 



On Fig. [71 we show for the two cases W = h and Vt^ = 80 the linear behavior in L of the averaged number of 
eigenstates < Ml > and of the variance Var{J\fL). Our final results concerning the compressibility as a function of 
the disorder strength W are 

W = 5,10,20,40,80 
X{W) = 0.45,0.45,0.62,0.79,0.87 (98) 

We find that our results for W = b and = 10 are nearly the same, as already found for the singularity spectrum 
/(a). Then the compressibility x grows with W as expected, and goes to 1 in the strong disorder limit W — > +00. 
A direct relation x + Di/d = 1 between the compressibility x of energy levels, and the information dimension 
Di = Uq^i of eigenfunctions has been recently conjectured and checked in various models [isj : for the present model, 
our numerical results do not seem compatible with this relation. For instance at = 5, the measured information 
dimension Di = ai ~ 0.15 would correspond via the conjectured relation in d = 1 to x' = 1 ^ Di ~ 0.85, whereas 
we measure the compressibility x — 0.45. Our conclusion is thus that the relation x + Di/d — 1 conjectured in [Tsj 
probably needs some hypothesis that is not satisfied by the present model. 



B. Results for the Cauchy distribution 

For the Cauchy disorder, our final results for the compressibility as a function of the disorder strength W read 

W =1,5,10,40 

X{W) = 0.60,0.64,0.77,0.96 (99) 
Here for W — 40, the obtained compressibility is very near the Poisson value 1. 



VIII. ANOMALOUS WEAK-DISORDER LIMIT 



A. Numerical results within the lower band 



In the previous sections IVI CI and IVIIi we have shown numerical results for the multifractal spectrum and the 
compressibility at the center of the band = for finite disorder W. However in the weak-disorder region, the 
density of states tend to break into two bands around the two pure delta peaks of Eq. |46] that are separated here 
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t'y ^7ven ~ ^odd ~ "^^1 ~ a consequcnce, the choice to work around E = is not appropriate anymore at weak 

disorder, but one can instead work in one of the two sub-bands. We have chosen to study the statistical properties of 
the eigenvalues and eigenvectors of a fraction 1/8 of the states of the lower sub- band (after checking that the density 
of states was nearly constant in this region). 

For the Box distribution, we find that the compressibility x takes the values 

W =0.1,1 

X^'^'^iW) - 0.41,0.42 (100) 
For the Cauchy distribution, we actually find the same limiting value 

^Cauchyfyf^ = 0.01) ~ 0.41 (101) 

These numerical results indicate that the compressibility remains finite in the weak disorder regime W 0^, whereas 
in other models, it vanishes smoothly in the disorder strength to recover the 'Random Matrix' value x = 0. 

We have also analyzed the multifractal spectrum of eigenstates for the values of W given above : they keep a 'strong 
multifractality' character, with a typical value around atyp — 2 and a minimal value amin near zero. 

B. Discussion 

The numerical results obtained in a given sub-band of the pure model indicate that the 'weak-disorder' regime is very 
anomalous in the Dyson hierarchical model. Indeed in usual models characterized by a non-degenerate continuum of 
plane waves, the non-degenerate perturbation theory yields the universal first order Gaussian correction in q{q— 1) to 
the multifractal spectrum (see [ISjim for more details on this 'weak multifractality' regime), and the compressibility 
is perturbatively close to the 'Random Matrix' value x = 0. Here these generic results do not apply, because the 
pure model is extremely degenerate (Eq. I46p . so that one should diagonalize the perturbation in each extensively 
degenerate subspace of the pure model. In some sense, this means that the perturbation is never 'weak', since there 
is no energy scale associated to the pure model in a given delta-peak. 

IX. CONCLUSION 

In this paper, we have described how the Dyson hierarchical model for Anderson localization, containing non-random 
hierarchical hoppings and random on-site energies, can reach an Anderson localization critical point presenting mul- 
tifractal eigenfunctions and intermediate spectral statistics, provided one introduces alternating signs in the hoppings 
along the hierarc hy ( instead of choosing all hoppings of the same sign as had been done up to now in the mathemat- 
ical literature Il3t-|l6ll ). This model is somewhat simpler than the 'ultrametric random matrices ensemble' considered 
by physicists [l7Hl9j, because here it is directly the matrix of non-random hoppings in each sample that presents 
a hierarchical block structure, whereas in References [l7l - [l9l |. it is only the matrix of the variances of the random 
hoppings that presents a hierarchical block structure. In particular, we have obtained exact renormalization equa- 
tions for some observables, like the renormalized on-site energies or the renormalized couplings to exterior wires. 
For the renormalized on-site energies, we have showed that the Cauchy distributions are exact fixed points. From 
the renormalized couplings to exterior wires, we have obtained that the typical exponent of eigenfunctions is always 
cttyp = 2 independently of the disorder strength, in agreement with our numerical exact diagonalization results for 
the box distribution and for the Cauchy distribution of the random on-site energies. We have also explained how 
this model has the same universal 'strong multifractality' regime in the limit of infinite disorder strength W — t- +oo 
as in other models. The big difference with other models is however that the singularity spectrum /(a) keeps for 
finite disorder a 'strong multifractal' character with very inhomogenous eigenfunctions, instead of flowing towards 
a 'weak multifractality regime'. The absence of this 'weak multifractality regime' comes from the anomalous pure 
spectrum of this hierarchical tree structure, with two extensively degenerate delta peaks (instead of some continuum 
corresponding to plane waves). 

We hope that the present work will stimulate further work in the renormalization analysis of this critical model, 
or in the 'critical ultrametric ensemble' which has only be studied via perturbation theory or numerics up to now 
since the main motivation to introduce Dyson hierarchical models is usually to obtain exact renormalization 
equations. A further goal is to better understand multifractality via renormalization in other non-hierarchical models. 
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as already discussed in Refs [H, HI] . 
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